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Abstract. A calculation of the current-quark mass dependence of hadron 
masses can help in using observational data to place constraints on the 
variation of nature's fundamental parameters. A hadron's u-term is a mea- 
sure of this dependence. The connection between a hadron's cr-term and the 
Feynman-Hellmann theorem is illustrated with an explicit calculation for the 
pion using a rainbow-ladder truncation of the Dyson-Schwinger equations: in 
the vicinity of the chiral limit a-,^ = 171^^/2. This truncation also provides a 
decent estimate of Up because the two dominant self-energy corrections to the 
p-meson's mass largely cancel in their contribution to a p. The truncation is 
less accurate for the w, however, because there is little to compete with an 
uj ^ pTT self-energy contribution that magnifies the value of fT(^ by < 25%. A 
Poincare covariant Faddeev equation, which describes baryons as composites 
of confined-quarks and -nonpointlike-diquarks, is solved to obtain the current- 
quark mass dependence of the masses of the nucleon and A, and thereby 
and a A- This "quark-core" piece is augmented by the "pion cloud" contribu- 
tion, which is positive. The analysis yields — 60MeV and a a ~ 50MeV. 



1 Introduction 

Chiral symmetry is explicitly broken in QCD by the current-quark mass term, 
which for the u- and d-quark sector is expressed in the action as 



where: {T^)ij = Sij and {r'^; k = 1,2, 3} are Pauli matrices; and m = (rriu+md) /2 
and rh = (m^ — m^l2. Empirical success with the application of chiral effective 
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theories to low-energy phenomena in QCD indicates that this term can often be 
treated as a perturbation. That simphfication owes fundamentally to the phe- 
nomenon of dynamical chiral symmetry breaking (DCSB) in QCD; namely, the 
feature that the dressed-quark Schwinger function is nonperturbatively modified 
at infrared momenta: p ^ 1 GeV. This is a longstanding prediction of Dyson- 
Schwinger equation (DSE) studies |^ HI 01 IH IS] that has recently been verified in 
numerical simulations of lattice-regularised QCD 0. A quantitative comparison 
and feedback between DSE and lattice studies is currently proving fruitful; e.g., 

Refs. [niHiiniiiniiniiiaiiiiiiiiisi. 

The (T-term for a hadron H is obtained from the isoscalar matrix element 

{H{x)\m Mz)\H{y)) , M^) = Qi^yQiz) , (3) 

and simple counting of field dimensions entails that for mesons the scalar form 
factor associated with this matrix element has mass-dimension two, while for 
fermions it has mass-dimension one. For all hadrons the u-term vanishes in the 
chiral limit. It is thus a keen probe of the impact of explicit chiral symmetry 
breaking on a hadron, in particular, as will be made plain below, on a hadron's 
mass. 

Such information is important for numerous reasons, some of longstanding 
jl6j . but our interest is prompted by its connection with the variation of nature's 
fundamental parameters |17| . It is a feature anticipated of models for the 

unification of all interactions that the so-called fundamental "constants" actu- 
ally exhibit spatial and temporal variation. In consequence there is an expanding 
search for this variation via laboratory, astronomical and geochemical measure- 
ments. An interpretation of some of these measurements can benefit from calcu- 
lations of the current-quark mass dependence of the parameters characterising 
nuclear systems. 

Of relevance herein are meson and baryon masses. A variation in light-meson 
masses will modify the internucleon potential, and a variation in the nucleon mass 
will affect the kinetic energy term in the nuclear Hamiltonian. Such changes could 
modify the binding energy in deuterium. That would have a material impact on 
Big Bang Nucleosynthesis (BBN) because the first step in BBN is the process 
p + n d + 7. The rate of this process is crucially dependent on the binding 
energy of deuterium, and this reaction is the seed for all subsequent processes 
and therefore the primordial abundance of light elements. A calculation of the 
current-quark mass dependence of hadron properties is therefore necessary to 
enable the use of observational data to place constraints on the variation of 
nature's constants |2nj . 

A variation of meson and nucleon masses will also modify the position of com- 
pound resonances in heavy nuclei {21^ . The position two-billion years ago of such 
a resonance in neutron capture by ^^^Sm has been determined using data from 
the Oklo natural nuclear reactor, with a non-zero shift reported cf. the present 
day [221 ■ The position of this resonance is very sensitive to the masses of mesons 
and the nucleon, and thus the Oklo data provide an acute method by which to 
measure the temporal rate of variation of nature's basic parameters. Data on the 
position of nuclear resonances several billion years ago are also available from the 
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study of nuclear reactions in stars Moreover, the interpretation of numer- 
ous measurements of quasar absorption spectra and superprecise atomic clocks in 
terms of the variation of nature's fundamental parameters also requires calcula- 
tions of the quark- mass-dependence of meson and nucleon properties |211 12311^ . 
These observations emphasise the important role that the calculations reported 
herein can currently play in the interpretation of many measurements performed 
in several areas of physics and astronomy. In addition, a large number of new 
and more accurate measurements are soon expected to appear. 

In Sect. 121 we focus on the cr-term for mesons, in particular, the vr, p and w. 
We present an explicit calculation of a-,^ from the pion's scalar form factor in 
the rainbow-ladder truncation of the DSEs. This is a useful means by which to 
exemplify the Feynman-Hellmann theorem in the context of QCD. In order to 
compute the c-terms of the nucleon and A it is first necessary to have a tool with 
which to calculate the masses of these baryons. Hence, in Sect.|31we recapitulate 
on a Poincare covariant Faddeev equation that serves this purpose In Sect.|l] 
we describe how this equation is used to calculate the dressed-quark contribution 
to (Tat and cr^, and how that result may be affected by meson-cloud corrections. 
Section [3 is an epilogue. It contains a tabulation of our "best estimates" for 
5mH/rnH in terms of 5rh/'fh. 



2 Mesons 



We judge that it is useful to begin with an explicit calculation of 0",^ from the 
pion's scalar form factor, which is defined by a momentum space re-expression 
of Eq. dSl): 

s.{Q'') = {^{P')\m MQ)\7r{P)) , = {P' - P)^ . (4) 

This three-point function arises, for example, in the analysis of vrvr scattering, as 
explicated in Ref. |26] . 

In the rainbow-ladder truncation of QCD's DSEs, 

sAQ^) = ircDF j-^S{l_^i)mr,,{l^i^o-,Q)S{l_^_i) 

2 ,v 2 '2 

where the trace is over colour, flavour and spinor indices, and 
^a,f3 = ^ + otP + f3Q. The rainbow-ladder approximation is the first term 
in a nonperturbative, systematic and symmetry preserving DSE truncation 
scheme \n\ I28j . This leading-order contribution preserves the one- loop renor- 
malisation group properties of QCD, and has provided a uniformly accurate 
description and prediction of a wide range of meson properties [2£:. (NB. We 
employ the Euclidean metric described in Appendix A: Euclidean Conventions.) 

Equation © is fully renormalised. It contains the renormalised rainbow- 
dressed quark propagator, which is the solution of 

S{p)-^ = Z^ii^-p + M'^n + ^ip), (6) 

/A y^a y^a 

9^Df,^{p - y7m'5(9) yTi^, (7) 



4 Sigma Terms of Light- Quark Hadrons 

wherein: represents a Poincare invariant regularisation of the integral, with A 
the regularisation mass-scale; D^^^ describes the interaction between light-quarks; 
and A4^^ describes the yl-dependent current-quark bare masses. The final step 
in any calculation is to take the limit A ^ oo. The quark-gluon- vertex and 
quark wave function renormalisation constants, Zi^2(C)^)) depend on the gauge 
parameter, the renormalisation point, and the regularisation mass-scale, but 
they are independent of the current-quark mass. Equation © is supplemented 
by the renormalisation condition 

S{p)-\2=^2=i7-P + M{0, (8) 

where A^(C) describes the renormalised (running) current-quark masses: 

Z^iC, A) M'^^'iA) = Z,{C, A)M{C) , (9) 

with Z4((^,yl) the Lagrangian-mass renormalisation constant. 

Equation © contains in addition: the Bethe-Salpeter amplitude for the 
ladder-pion bound state, which is obtained from 

r^k; P) = - j\^D^,{k - q)^^^S{q+))r^{q; P)S{q^)^^,, (10) 

with q± = q^ -P/2, and normalised canonically 

= ticDF jy.{q]-P)^S{q + P/2)r^{q-P)S{q-P/2) 

+ ^'^CDF r^{q--P)S{q + P/2)r^{q-,P)^S{q-P/2)- (11) 

and the isoscalar vertex 

r,o(£;g) = Z4(C,^)r°+^V^/..(^-'?)Y7/.5((/+)r,o(£;Q)5(g_)y7.- (12) 

NB. Although the isoscalar vertex r^o{l;Q) depends on the renormalisation 
point, the product fh{Q r^o{l;Q) is renormalisation-point-independent. 
The pion cr-term is defined via 

2m^a^ := s^(g2 = q) (13) 
r d'^f 

= ticDF J 5(^_i,o) fh r,o (£_i,o; 0) 5(^_i,o) 
Equation H14|l can be simplified using a textbook result; viz.. 



S{k) = -S{k)r^o{k;0)S{k), 



(15) 
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and hence 



2 a. = - m(c) trccF / ■ ,„; -P) £(() P) 



dm{C) dm{() 
where the last hne fohows from Eq. so that 

We emphasise that in arriving at Eq. (|16() we have depended heavily upon 
the fact that the rainbow-ladder expression is the leading term in a systematic, 
nonperturbative and symmetry preserving truncation of the DSEs |1U1 1271 05] . 
The derivation provides a concrete illustration of a general result that may be 
viewed as a consequence of the Feynman-Hellmann theorem. Applied in the 
present case, this theorem states that the response of an eigenvalue of the QCD 
mass-squared operator to a change in a parameter in that operator is given by 
the expectation value of the derivative of the mass-squared operator with respect 
to the parameter. The result is valid in this form for all mesons; i.e., 

^muOM ■■= sm{0) = "^iO-Q^^ ^ = "^(0^^=^ • (17) 

NB. The cr-term is a renormalisation point invariant, in general and also in our 
explicit calculation. 

In QCD the pion's mass is expressed precisely via pUj : 

ml = -2m{0^- (18) 
In this expression, is the pion's leptonic decay constant: 

U5'^P^ = Z2iC,A)tr [^yro-/^S{q+)niq;P)S{q^); (19) 

Jq 

and is the residue of the pion pole in the dressed-quark-antiquark-pseudoscalar 
vertex:^ 

ip^OS"'=Z,{C,A)tvl yj,S{q+)rUq;P)S{q-). (20) 



^NB. The factor of Z2{C,A) in Eo. 11911 guarantees that the right-hand-side (r.h.s.) is gauge 
invariant, independent of the renormalisation point and finite in the hmit A oo; and ^4(C, A) 
in Eq. I|2()^ ensures that the r.h.s. is gauge invariant, finite when the regularisation scale is 
removed, and evolves according to the renormalisation group precisely as required to guarantee 
that the product rh{(^)pJ^Q is renormalisation point independent. 
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Additional information may now be obtained from Eq. by substituting 
Eq. H18|l . In the neighbourhood of the chiral Hmit ^30j 

PAC) = jQ-, (21) 

where (qq)^ is the vacuum quark condensate and is the chiral limit value of 
the leptonic decay constant, and hence 

„ m-O „ _ _^ m~0 1 

2m^CT^ = -^m{C) =^ (Tn = -m^- (22) 

(/tt) ^ 

This is an essential consequence of DCSB. 

With typical, calculated chiral limit values of [5^: = 0.088 GeV and 
{QQ)(=iGeY ~ (— 0.241 GeV)^, and a current-quark mass m(C = 1 GeV) = 
0.0055 GeV, one estimates from Eq. 

= 0.95 ^/m{C = 1 GeV) = 0.071 GeV . (23) 

On the other hand, the value obtained via direct calculation in rainbow-ladder 
truncation can be determined from the information provided in Refs. |311l32ll^ 
using Eq. (|T6|) : namely, 

(j^^ = 0.069 GeV. (24) 

It is useful and relevant to compare the pion cr-term with that obtained 
for the p-meson, and one may infer the rainbow-ladder truncation result from 
Refs. ESI; viz., 

2 mp af^ = s^^(O) = (0.194 GeV)^ ^ af^ = 0.025 GeV . (25) 

Since the p-meson is unstable it is natural to ask whether meson-loop self- 
energy corrections alter this result significantly. It was shown in Ref. |S1] that the 
p — > TTTT and p UJTT self-energy insertions, 11 and 11 ^'^ respectively, provide 
the largest shift in the mass of the p-meson. Hence, we focus on the contribu- 
tion of these terms to dp and estimate their effect by adapting the approach of 
Ref. |SSI) which gives formulae for Up'" and 11 ^'^ that can be used to calculate 
their m^-dependence. For all values of ttt-tt, both 11^'^ < and U^^ < 0; viz., 
they both act to reduce m^. However, since the p — > tttt decay channel is open 
in the vicinity of the experimental value of uIt^, we find 



0<ml ^ 



dml 



QJJTTTT 

= (0.056GeV)2<-m2— ^ 



= (0.090 GeV)^ 

(m|)expt. 

(26) 

Hence the contributions to cjp from these self-energy insertions cancel each other 
to a large degree, and one obtains 

2 nip (7^^+^^+'^^ = (0.181 GeYf a^^+^^^^^ = 0.022 GeV . (27) 

This suggests that in connection with Up the effect of meson-loop contributions is 
moderate; i.e., < 15 %, so that the rainbow-ladder truncation yields a fairly good 
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estimate. It is thus noteworthy that with the interaction employed in Ref. |S21 
the rainbow-ladder truncation yields, via the formulae presented in Ref. [HH] . 

cj^^ = 0.034 m^-^ . (28) 

At this point we have sufficient information to estimate the w-meson's a- 
term. In the case of isospin symmetry, the uj- and p-mesons are indistinguishable 
in rainbow-ladder truncation .^(i and hence 

a^-^ = a^^. (29) 

However, the states differ in the nature of their self-energy corrections. In contrast 
with the /3-meson, the u — > vrvr coupling vanishes. This leaves only the to ^ pir 
self-energy, and in the isosymmetric limit^ 

nP]" = n'^'' . (30) 

From this and Eq. (|26() 

2 a^^+P^ = (0.217 GeV)^ a^^+P^ = 0.032 GeV = 0.043 , (31) 

which is an increase of < 25% over the rainbow-ladder result. The marked dif- 
ference between the values in Eqs. (|28)) and (|31|) . and that in Eq. (|22j) highlights 
again the particular character of QCD's Goldstone mode. 



3 Baryons 

In a consideration of Eq. Q and the cr-term for fermions, the analogue of Eq. ()17|) 
is 

ap:=sH0)=m{O^y (32) 

This is apparent in the context of the Feynman-Hellmann theorem because 
fermion fields have mass-dimension 3/2 and their spectrum is described by a 
mass operator, not a mass-squared operator. The cr-terms of the nucleon and 
A can thus be determined in any framework that provides for a calculation of 
the current-quark-mass dependence of Mjy and M^. We will use the Poincare 
covariant Faddeev equation model described in Refs. |i25j . 



3.1 Faddeev Equation 

In order to explain the input it is useful to recapitulate on the model's foundation. 
For quarks in the fundamental representation of colour-S'[/(3): 

3c 3c 3c = (3c e 6c) 3c = Ic e 8'c 8c © 10c , (33) 

and hence any two quarks in a colour-singlet three-quark bound state must con- 
stitute a relative colour- antitriplet. This fact enables the derivation of a Faddeev 

^NB. The LO — > tvtt-k contribution is negligible ^37^ and no other strong decay channels are open. 
The signs are incorrect in Eqs. (2.11) and (3.4) of Ref. |37| . 



8 



Sigma Terms of Light- Quark Hadrons 



equation for the bound state contribution to the three quark scattering kernel 
PH] because the same kernel that describes mesons so well [2^1 is also attractive 
for quark-quark scattering in the colour-3 channel. 

In this truncation of the three-body problem the interactions between two 
selected quarks are added to yield a quark-quark scattering matrix, which is then 
approximated as a sum over all possible diquark pseudoparticle terms |39| : Dirac- 
scalar + -axial- vector +[...]. The Faddeev equation thus obtained describes the 
baryon as a composite of a dressed-quark and nonpointlike diquark with an 
iterated exchange of roles between the bystander and diquark-participant quarks. 
The baryon is consequently represented by a Faddeev amplitude: 

<F = + ^2 + '^'3 , (34) 

where the subscript identifies the bystander quark and, e.g., are obtained 
from !P'3 by a correlated, cyclic permutation of all the quark labels. 

3.2 Ansdtze for the Nucleon and A 

We employ the simplest realistic representation of the Faddeev amplitudes for 
the nucleon and A. The spin- and isospin-1/2 nucleon is a sum of scalar and 
axial-vector diquark correlations: 

^3ip^,ai,n)=^^t (35) 

with (pi, ai,Ti) the momentum, spin and isospin labels of the quarks constituting 
the bound state, and P = P1+P2+P3 the system's total momentum.^ Since it is 
not possible to combine an isospin-0 diquark with an isospin-1/2 quark to obtain 
isospin-3/2, the spin- and isospin-3/2 A contains only an axial- vector diquark 
component 

W^{pi,ai,n)=Vl''. (36) 
The scalar diquark piece in Eq. (|35j) is 

Mfip^,a,,n) = [r'\lp[,2Y,K)]l\ZA'^iK)[Sii-,P)n{P)]2\, (37) 

where: the spinor satisfies (recall [Appendix A:| Euclidean Conventions) 

(^7 • P + M) u{P) = = u{P) {i-f-P + M), (38) 

with M the mass obtained by solving the Faddeev equation, and it is also a 
spinor in isospin space with ip^ = col(l,0) for the proton and ip^ = col(0, 1) 
for the neutron; K = pi + p2 =: P{i2}, P[i2] = Pi - P2, ^ ■= (-^{12} + 2p3)/3; 
A^^ is a pseudoparticle propagator for the scalar diquark formed from quarks 
1 and 2, and P'^^ is a Bethe-Salpeter-like amplitude describing their relative 
momentum correlation; and S, a 4 x 4 Dirac matrix, describes the relative quark- 
diquark momentum correlation. (5, P^^ and A^^ are discussed in Sect. 13.31 ) The 

^NB. Hereafter we assume isospin symmetry of the strong interaction; i.e., the u- and d-quarks 
are indistinguishable but for their electric charge. This simplifies the form of the Faddeev 
amplitudes. 
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colour antisymmetry of l^'a is implicit in F , with the Levi-Civita tensor, ecic2C3, 
expressed via the antisymmetric Gell-Mann matrices; viz., defining 

{H^ = iX\ = -iX\ = iX^} , (39) 

then €cic2C3 = {H^')c,c2- [See Eqs.®, (EZI)-] 
The axial-vector component in Eq. (|35() is 

Af'\p,,a,,n) = [t^rl,\lp[,,y,K)]ZZKtiK)[Al{e-,P)u{P)]Z, (40) 
where the symmetric isospin-triplet matrices are 

and the other elements in Eq. H4U|) are straightforward generalisations of those 
in Eq. (jSZI). 

The general form of the Faddeev amplitude for the spin- and isospin-3/2 A is 
complicated. However, isospin symmetry means one can focus on the A~^^ with 
it's simple flavour structure, because all the charge states are degenerate, and 
consider 

Pr = [t+rf (ip[i2];i^)];\^44:(K)[P.p(^;P)np(P)^+];='3, (42) 

where Up{P) is a Rarita-Schwinger spinor, Eq. (|A.10|) . 

The general forms of the matrices S{i;P), A1{i; P) and 'D^p{£;P), which 
describe the momentum space correlation between the quark and diquark in the 
nucleon and the A, respectively, are described in Ref. j4Uj . The requirement that 
S{i;P) represent a positive energy nucleon; namely, that it be an eigenfunction 
of A+{P), Eq. (|A.7|) . entails 

S{i;P)=si{e;P)lB + (il-i-i-PlD) S2{i;P), (43) 

where (/D)rs = Srs, (-"^ = 1, P^ = —1. In the nucleon rest frame, si^2 describe, 
respectively, the upper, lower component of the bound-state nucleon's spinor. 
Placing the same constraint on the axial- vector component, one has 

6 

P)=Y1 Pni^' ^) ^5 ^"(^; P) , i = +, 0, - , (44) 
n=l 

where (£> = 4 + | • P P^, 7^ = 7^ + 7 • P P^) 



Al = j-i^P,, Al = -iP,, Al = j-i^i\ 

At = J A^ — Iv ~ Af, , Af, = i^^, 7 • £ — A^ . 



(45) 



Finally, requiring also that T)yp{£] P) be an eigenfunction of yl+(P), one obtains 
V,p{l- P) = S^ii; P) + 75^^(^; P) if , (46) 
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with S'^ and given by obvious analogues of Eqs. and respectively. 
One can now write the Faddeev equation satisfied by 'F^ as 



S{k;P)u{P) 
A^^{k;P)u{P) 



The kernel in Eq. H47|) is 

M{k,i; P) 

with 



(2vr)' 



■Mik,i;P) 



S{t,P)u{P) 
Ai{t,P)u{P) 



(47) 



(48) 



M 



00 



{kq ^qq/'^]^qq)S {iqq nq 



kq) i^q kqq/2; kqq) S{lq) ^ i^qq) > 

(49) 

where: £q = £ + P/3, kq = k + P/3, Iqq = -l + 2P/3, kqq = -k + 2P/3 and the 
superscript "T" denotes matrix transpose; and 



{M,l)i=t^rl\kq-£qq/2-lqq) 



' (^gg kq) r 
{Ml0)'=r'\kq-£qq/2;eqq) 



{iqq kq) {£q kqq/2; kqq) S{£q) Apj^{£qq) , 



xS'^{£qq kq)t^r^ (Iq kqq/2; kqq) S{£q) i^qq) i 



{Ml^)%=t^rl\kq-£qq/2;£qq) 
S {£qq ~ k 

The Z\'s Faddeev equation is 



xS'^{£qq kq)t^r^ {£q kqq/2; kqq) S {£q) A^j^ (£qq) . 



(50) 
(51) 
(52) 



V,,{k;P)u,{P) = aJ ^Mt^{k,£;P)V,^{£;P)u^{P), (53) 



with 



J^\^ — t^-f^CT {kq — ('qq/'^',iqq) S^{£qq — kq)t^rl {£q — kqq/2; —kqq) S{£q) A\^{£qq) . 

(54) 



3.3 Completing the Faddeev Equation Kernels 

To complete the Faddeev equations, Eas. (|17l) & one must specify the 

dressed-quark propagator, the diquark Bethe-Salpeter amplitudes and the di- 
quark propagators that appear in the kernels. 



3.3.1 Dressed-quark propagator 

The dressed-quark propagator has the general form 

S{p) = -ij.pav{p'^)+as{p'^) = l/[i-fpA{p'^)+B{p'^)] 



(55) 
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and can be obtained from QCD's gap equation, the rainbow-ladder truncation of 
which is given in Eq. It is a longstanding prediction of DSE studies in QCD 
that the wave function renormalisation and dressed-quark mass: 

= ^ M{p^) = , (56) 

respectively, receive strong momentum-dependent corrections at infrared mo- 
menta j^miSl: Z{p'^) is suppressed and M(p'^) enhanced. The enhancement of 
M{p^) is central to the appearance of a constituent-quark mass-scale and an 
existential prerequisite for Goldstone modes. The mass function evolves with 
increasing to reproduce the asymptotic behaviour familiar from perturbative 
analyses, and that behaviour is unambiguously evident for > lOGeV^ [I H II ^ . 

The impact of this infrared dressing on hadron phenomena has long been 
emphasised |4,S| and, while numerical solutions of the quark DSE are now readily 
obtained, the utility of an algebraic form for S{p) when calculations require the 
evaluation of numerous multidimensional integrals is self-evident. An efficacious 
parametrisation of S{p), which exhibits the features described above, has been 
used extensively in hadron studies El EH] • It is expressed via 

as{x) = 2mJ^{2{x + rh'^)) + J^{bix)J^{bsx) [bo + b2J^{ex)] , (57) 
av{x) = [1 - .F(2(x + m^))] , (58) 

with x = p^ /}?, fh = m/A, 

1 - e~^ 

T{x) = , (59) 

X 

as{x) = XcFsip^) and ay(x) = A^cjy(p^). The mass-scale, A = 0.566 GeV, and 
parameter values^ 



m bo bi 62 bs 
0.00897 0.131 2.90 0.603 0.185 



(60) 



were fixed in a least-squares fit to light-meson observables [441 145j . The dimen- 
sionless u = d current-quark mass in Eq. (|En|) corresponds to 

m = 5.08MeV. (61) 

The parametrisation yields a Euclidean constituent-quark mass 

M^^ = 0.33 GeV, (62) 

defined as the solution of = M^(p^) |31j. Since the dressed-quark mass function 
is a renormalisation point invariant then so is M^. Hence, one may define a 
constituent-quark cr-term 

^Q:="^(C)^ = 6.2MeV (63) 



*e = 10 in Eq. 1571 acts only to decouple the large- and intermediate-p^ domains. 
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for the parametrisation we have just described. NB. With the model interaction 
of Ref. j32j . Eq. (jH)) gives ctq = 9.4 MeV in rainbow-ladder truncation; and that 
of Ref. ^01 gives ctq ~ 8- 9 MeV, depending on the Ansatz for the quark-gluon 
vertex, with the lower value corresponding to rainbow-ladder truncation. It is 
thus likely that the parametrisation of S{p) employed herein underestimates aq, 
owing to the implicit assumption that &o,i,2,3 Stice m-independent. 

3.3.2 Diquark Bethe-Salpeter amplitudes 

The rainbow-ladder DSE truncation yields asymptotic diquark states in the 
strong interaction spectrum. Such states are not observed and their appearance 
is an artefact of the truncation. Higher order terms in the quark-quark scattering 
kernel, whose analogue in the quark-antiquark channel do not much affect the 
properties of vector and flavour non-singlet pseudoscalar mesons, ensure that 
QCD's quark-quark scattering matrix does not exhibit singularities which corre- 
spond to asymptotic diquark states jlOj. Nevertheless, studies with kernels that 
do not produce diquark bound states, do support a physical interpretation of 
the masses, ^{qq)jp , obtained using the rainbow-ladder truncation: the quan- 
tity l(qq)jp = l/m(gg)^^ may be interpreted as a range over which the diquark 
correlation can persist inside a baryon. These observations motivate the Ansatz 
for the quark-quark scattering matrix that is employed in deriving the Faddeev 
equation: 

[Mqqik, q; = Yl ^''"(^^ ^^"(^) ^■'"(9; K) . (64) 

jp=o+,i+,... 

One practical means of specifying the F"^^ in Eq. H64|l is to employ the so- 
lutions of a rainbow-ladder quark-quark Bethe-Salpeter equation (BSE). Using 
the properties of the Gell-Mann matrices one finds easily that := F'^ C'^ 
satisfies exactly the same equation as the colour-singlet meson hut for a 
halving of the coupling strength jlS]. This makes clear that the interaction in 
the 3c {qq) channel is strong and attractive. For the correlations relevant herein, 
models typically give masses (in GeV) |^ HH] : 

^{ud)^+ = 0.74 - 0.82 , m(„„)^^ = m(„d)^_^ = m(^dd)^+ = 0.95 - 1.02 . (65) 

Such values are confirmed by results obtained in simulations of quenched lattice- 

QCD gg. 

A solution of the BSE equation requires a simultaneous solution of the quark- 
DSE j48] . However, since we have already chosen to simplify the calculations by 
parametrising 5(p), we also employ that expedient with F'^ , using the following 
one-parameter forms: 

F^\k-K) = -L^H'^Ct^^iT^Hk^/u^l^), (66) 

t'Fl\k-K) = -}^H^^J^Ct^J^{k'/u;j^), (67) 

with the normalisation, Af'^'', fixed by an appropriate analogue of Eq. I|1H) |25| . 
These Ansdtze retain only that single Dirac-amplitude which would represent a 
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point particle with the given quantum numbers in a local Lagrangian density: 
they are usually the dominant amplitudes in a solution of the rainbow-ladder 
BSE for the lowest mass diquarks [HlllS] and mesons |3I1 EH • 

3.3.3 Diquark propagators 

Solving for the quark-quark scattering matrix using the rainbow-ladder trun- 
cation yields free particle propagators for A"^ in Eq. (|64() . As already noted, 
however, higher order contributions remedy that defect, eliminating asymptotic 
diquark states from the spectrum. The attendant modification of Z\ can be 
modelled efficiently by simple functions that are free-particle-like at spacelike 
momenta but pole-free on the timelike axis ^U]; namely,^ 

A'\k) = -^^{K'/u^l^), (68) 

<W = U. + ^)^HKV^I.), (69) 

where the two parameters rrijp are diquark pseudoparticle masses and ojjp are 
widths characterising . Herein we require additionally that 

= 1 ^ u;2^ = 1 , (70) 

ii'2=0 

which is a normalisation that accentuates the free-particle-like propagation char- 
acteristics of the diquarks within the hadron. 

3.4 Nucleoli and A Masses 

All elements of the Faddeev equations, Eas. (|17|l & ((^ . are now completely 
specified. We solve the equations via the method described in Ref. [H^. Owing to 
Eq. H70() , the masses of the scalar and axial- vector diquarks are the only variable 
parameters. The axial- vector mass is chosen so as to obtain a desired mass for 
the A, and the scalar mass is subsequently set by requiring a particular nucleon 
mass. 

Two primary parameter sets are presented in Table ^ Set A is obtained 
by requiring a precise fit to the experimental nucleon and A masses. It has long 
been known that this is possible; e.g., Ref. reports octet and decuplet baryon 
masses in which the rms deviation between the calculated mass and experiment is 
only 2%. However, it is also known that such an outcome is undesirable because, 
e.g., studies using the cloudy bag model indicate that the nucleon's mass 
is reduced by as much as SMn = —300 to — 400MeV through pion self-energy 
corrections |57| . Furthermore, a perturbative study, using the Faddeev equation, 
of the mass shift induced by pion exchange between the quark and diquark 

^These forms satisfy a sufficient condition for confinement because of the associated violation 
of reflection positivity. This notion may be traced from Refs. |51l 1521 l53l 154) and is reviewed 
in Refs. ElllE]. 
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Table 1. Mass-scale parameters (in GeV) for the scalar and axial- vector diquark correlations, 

fixed by fitting nucleon and A masses: Set A provides a fit to the actual masses; whereas 

Set B provides masses that are offset to allow for "pion cloud" contributions (59) . We also 

list ujjp = ——nijp, which is the width-parameter in the {qq)jp Bethe-Salpeter amplitude, 
v2 

Eqs. fc ll^: its inverse is an indication of the diquark's matter radius. Sets A* and B* 
illustrate effects of omitting the axial-vector diquark correlation: the A cannot be formed and 
Mjv is significantly increased. It is thus plain that the axial- vector diquark provides significant 
attraction in the Faddeev equation's kernel. 



set 


Mn 


Ma 


mo+ 


mi+ 






A 


0.94 


1.23 


0.63 


0.84 


0.44=1/(0.45 fm) 


0.59=1/(0.33 fm) 


B 


1.18 


1.33 


0.79 


0.89 


0.56=1/(0.35 fm) 


0.63=1/(0.31 fm) 


A* 


1.15 




0.63 




0.44=1/(0.45 fm) 




B* 


1.46 




0.79 




0.56=1/(0.35 fm) 





constituents of the nucleon obtains 6Mn = —150 to — 300MeV (SS]. This leads 
to Set B, which was obtained by fitting to nucleon and A masses that are inflated 
so as to allow for the additional attractive contribution from the pion cloud |59j . 

4 Nucleon and A fj-terms 

4.1 Analysis 

In order to calculate a^^^ in the framework we have outlined it is necessary 
to know the variation with current-quark mass of the elements described in 
Sects. l3?3?Tt43.3.3l That is straightforward for the dressed-quark described in 
Sect. l3.3.H so long as we assume that the explicit m-dependence is dominant in 
the neighbourhood of the physical current-quark mass. As we saw in connection 
with Eq. (|63|) . this assumption may underestimate the m-dependence. 

The diquark propagators in Sect. 13. 3^ each involve a single mass-scale, and 
the response of the nucleon and A masses to changes in these mass-scales is much 
as one would expect if niQ+ & mi+ were true diquark masses. This is apparent 
in Tabled For example, in changing from Set A to B, M^v increases by 240 MeV 
when (777-0+ + ^^1+) is raised by 210 MeV; and the increase of 1711+ by 50 MeV 
raises by 100 MeV. This is an appealing feature of the model. 

However, the connection should not be viewed as trivial because, in addition, 
these mass-scale parameters appear in the diquark Bethe-Salpeter amplitudes 
described in Sect. 13. 3^ These amplitudes constitute the final element in the 
Faddeev equation kernel. In principle their pointwise form will evolve with the 
current-quark mass and in our model, owing to Eq. (|7(H) . that is effected through 
an explicit dependence on the mass-scale characterising the appropriate diquark's 
propagator. This entails that the momentum-space width grows with increasing 
diquark mass-scale, a property which is consistent with the Bethe-Salpeter equa- 
tion studies of Ref. [HJ- It follows that the normalisation factor in Eqs. (|66|) and 
1)6 7|) becomes larger with increasing diquark mass. That leads to a reduction 
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in the quark-quark <-> diquark transition probability amplitude and thereby a 
decrease in the amount of attraction in the Faddeev equation's kernel. This ex- 
plains why, in changing from Set A to B, and M/^ rise by more than just the 
increase in the diquark masses. 

Diquark masses can unambiguously be defined and calculated in rainbow- 
ladder truncation, and that calculation would provide the variation of mQ+ and 
mi+ with current-quark mass; viz., an estimate of the diquark fi-terms (To+ and 
ai+. Although a realistic calculation of mQ+ and mi+ is available the evolu- 
tion of these masses with current-quark mass is not. We have therefore studied 
this problem using the algebraic model of Ref . j6Uj , which gives (in MeV) 

ao+ = 21>ai+ =20>ap = 17; (71) 

i.e., an evolution with current-quark mass of a diquark's mass which is similar 
to that of the /j-meson. This being the case, we proceed by supposing that the 
current-quark mass dependence of mo+ and mi+ is well approximated by that 
of the /3-meson; namely, we adapt the fit of Ref. jHH] as follows: 

mo+(m) = 771°+ -I- ^CQva + cim , (72) 
mi+{m) = + ^/c(ym + cim , (73) 

with Co = 0.195 GeV, ci = 1.90, and (in GeV) 





Set A 


Set B 




0.59 


0.75 




0.80 


0.85 



This assumption, too, will likely bias our calculation toward an underestimate 
of (TAT and aA- 

Now that the response of the Faddeev equation's kernel to changes in 
the current-quark mass has been elucidated, the calculation of and a a 
is straightforward. It is only necessary to evaluate M^r and M/^ on a small 
domain of current-quark mass, centred on the physical value, Eq. (jblf) : viz., 
4.5 < m (MeV) < 5.5, and interpolate. The a-terms follow from Eq. (|32|) and 
our results are presented in Table |2J It is notable that so long as Eq. ((7T|) is valid. 

With the kernel we have employed, the Faddeev amplitudes describe a 
baryon's dressed-quark core. Differences between the amplitudes obtained with 
Set A and Set B are exposed in calculations of the nucleons' electromagnetic 
form factors ^25 . It is apparent in Table [J] that the values of the cj-terms are 
almost independent of whether Set A or Set B is used. For the nucleon, the value 
is influenced much more by whether, or not, the axial-vector diquark is retained. 
It is thus evident that ajy and are primarily determined by the evolution 
of the constituents' masses with current-quark mass, and that the pointwise re- 
sponse of the Faddeev amplitude affects the results very little. In this connection 
we observe that the cr-term is the = value of a form factor and hence is 
determined by global rather than local properties of the distribution of bound 
state constituents. 
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Table 2. "Quark core" a-terms for nucleon and A evaluated from Faddeev equation solutions 
using Eq. 13211 . All dimensioned results tabulated are reported in GeV. NB. On the domain 
4.5 < m (MeV) < 5.5, Mjv(m) and AlA{m) are linear. The Set label indicates which parameters 
from Table were used in solving the Faddeev equation. The row labelled SetBg- is described 
in connection with Eq. I|82|l . Recall that absent an axial- vector diquark there is no A. For 
comparison, a value of (tat may be inferred from the isospin-even elastic irN scattering amplitude 
and analyses yield: an ~ 45 MeV i62||nS| or o-]v ~ 67 ± 8MeV ElinS]. Moreover, analyses of 
data from simulations of two-flavour lattice-regularised QCD yield: itjv ~ 51 MeV |52l; <^n ~ 
51 . . . 54 MeV gZl; and cr]v ~ 49 MeV and a a ~ 21 MeV _68.. A chiral-quark model analysis 
gives (jjv — 45 ±5 MeV and a a ~ 32 ±3 MeV ^69^. Our best estimates are presented in Eq. I|83|l . 



Set 


Mn 


Ma 




0"Z\ 


A 


0.94 


1.23 


0.046 


0.042 


B 


1.18 


1.33 


0.047 


0.042 


A* 


1.15 




0.040 




B* 


1.46 




0.041 






1.18 


1.33 


0.057 


0.051 



Here, too, it is natural to ask whether meson-loop self-energy contributions 
materially affect the results in Table|2l This may be addressed following Ref. j59j . 
which explains the nature of the vrA^-loop corrections to the Faddeev equation 
nucleon mass. A reliable estimate of this mass shift is 



5Mn = - 67r / — ' — ^ , (75) 



where /^^^ = g%j^^ml/{16TrMff), uj^ik"^) = ^PTrnf, uJNik'^) = ^^k^ + mlj, 
and, with a Pauli-Villars regularisation of the nucleon's pion-induced self-energy, 

tx(fc2) = i/(i + fc2/A2). (76) 

In Eq. ()76() . A is the pion-loop regularisation scale. It is finite and nonzero in a 
realistic calculation because neither the pion nor the nucleon is pointlike. 

To estimate the contribution to cr at from the mass shift in Eq. ((TSJ we evalu- 
ated this expression using a fixed value of {/nn-k/Mn) = {qa/ I-k) = 13.6 GeV~^, 
and values of Mn and m-7r that vary with the current-quark mass, m; namely, 
we used the nucleon mass obtained from the Set B Faddeev equation at a given 
value of m, which is described by (in GeV) 

= 1.13 + 9.20 m , m G [0.0045, 0.0055] (77) 

and [SOI 

m^ = -2m ^ (78) 

J TT 
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where, with the model of Sect. nmTl the calculated values are [IK] : 
{qqY = (-0.25 GeV)^ and = 0.090 GeV. This procedure gave (in GeV) 



A 


0.30 


0.40 




0.0037 


0.0058 



(79) 



viz., a contribution to from the vrA^ self-energy diagram of ~ 4 — 6MeV. NB. 
We list results with values of the regularisation parameter, A, that are typical of 
realistic analyses |251 161| . The smaller value of A is favoured in our model j25j . 
Moreover, the values in Eq. (|7n|l are unchanged if the Set A Faddeev equation 
results for Mjy are used. 

Comparing rows 2 & 4 and 3 & 5 in Table 121 it is apparent that the vr A'^-loop 
contribution to ctat, Eq. H79|). possess the same sign and is similar in magnitude 
to that arising from the axial-vector diquark. This is typical of the constructive 
interference between the pion cloud and axial- vector diquark correlations in the 
Faddeev equation. 

The shift in the nucleon mass owing to a ttA self-energy contribution has the 
same sign as that in Eq. ()75() and is no larger in magnitude. This suggests that 
the total correction to cjat from vrA^ and ttA self-energy corrections is 

~ 7MeV. (80) 

The shift in owing to vrA^ and ttA self-energy corrections is approximately 
one-half of the analogous correction to Mpf |70j . (Naturally, this fact lays behind 
the choice of SetB parameters in Tabled) Hence we deduce 

~3MeV. (81) 

The primary remaining source of uncertainty within our model is the assump- 
tion associated with Eqs. ((7^ . H73() . which entails o"o+ = cri+ = a p. To develop 
a notion of the error this may introduce we repeated the calculation of cttv and 
a A using the SetB Faddeev equation parameters with aQ+ = ai+ = 1.2 ap and 
the ^nsatz: 

mgg{m) = + m cr^g , (82) 
with rrigg given in Eq. (|74|) . This gave the results listed as Set B^ in Table IS 

4.2 Result 

Combining the elements discussed above we judge that, in the model employed, 
the best estimate of the cr-terms is 

o-AT ~ 60 MeV , cj^ ~ 50 MeV . (83) 

5 Epilogue 

We illustrated the connection between the cr-term and the Feynman-Hellmann 
theorem with an explicit calculation of cTtt in the rainbow-ladder truncation of 
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the Dyson-Schwinger equations. Naturally, in the vicinity of the chiral limit we 
obtain the model-independent result: cr^ = 

While on the subject of mesons, we also calculated dp, for which the rainbow- 
ladder truncation likewise provides a reliable estimate. A priori that need not 
have been the case because p tttt and p lott self-energy terms noticeably 
reduce m^. However, while the contribution to ap from each term is relatively 
large, they are of similar magnitude but opposite sign and therefore cancel to 
a large extent. We find dp ~ mp/34, a value that is appreciably smaller than 
Cat O E2l- For the w-meson, on the other hand, there is little to cancel the 
contribution to a^j from uj ^ pn, which is three-times larger than that to ap from 
p — > uJTT. Hence a^j — muj/23. A comparison of these values with cr,r highlights 
the importance of dynamical chiral symmetry breaking to light-hadron masses. 

To calculate and a a we employed a Poincare covariant Faddeev equation, 
which describes baryons as composites of confined-quarks and -diquarks. Two 
parameters appear in the model Faddeev equation. They are the mass parameters 
of the scalar and axial-vector diquark correlations, and were fixed by fitting 
stipulated masses of the baryons. The dependence of these diquark mass-scales 
on the current-quark mass was a key assumption in our calculation of the baryon 
cr-terms and the principal source of uncertainty. 

The Faddeev equation solution describes a baryon's "quark core" . This must 
be augmented in a consistent fashion by chiral-loop corrections. Those chiral 
corrections increase ajy by < 15% and by < 10%. 

Our analysis indicates that cjat is large compared with estimates based on 
chiral effective theory. Nonetheless, our value of aN ~ 0.06 Mjv is consistent 
with modern experimental evaluations j(i4| 165) . 

It has been argued, based on a perturbative treatment of SU (3)-flavour sym- 
metry breaking [231 ? that a large value of implies a large value for the ratio 
y = 2{p\ss\p) / {p\uu + dd\p); where large means y/2 > 0.2. While we have not 
studied this ratio, there are numerous examples of explicit calculations in which 
that inference is invalid because baryon masses exhibit a strong nonlinear de- 
pendence on the s-quark mass; e.g., Refs. (69l W^ 175 1 176]. 

With this study we have gained insight into the current-quark mass depen- 
dence of hadron properties. It is of use to gather our numerical results and 
reiterate them here in the form 

SruH (TH 6m 

— — , l,o4j 

where fh is defined in Eq. ((2)). At the physical value of the light current-quark 
mass we have 





IT p 


to 


N 


A 




0.498 0.030 


0.043 


0.064 


0.041 


rriH 











(85) 



From this point one may proceed directly to investigate the impact of our re- 
sults on the analysis of experimental observations that relate to the variation of 
nature's fundamental "constants." 
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Appendix A: Euclidean Conventions 

In our Euclidean formulation: 

4 

p - q = ^Piqi; (A.i) 

i=l 

{Tm. 71-} = '2Sf,^; 7m = 7/^ ; (^ni^ = ^ [7m> 7.'] ; tr[767M7i,7p7a] = -4 e^,^pa , ei234 = 1 - (A.2) 
A positive energy spinor satisfies 

u{P, s)ii-yP + M) = 0={i'yP + M) u{P, s) , (A.3) 
where s = ± is the spin label. It is normalised: 

u{P, s) u{P, s) = 2M (A.4) 

and may be expressed explicitly: 

u{P, s) = VM-i£ ( ^a P ] , (A.5) 
V M-i£^' J 

with £ = + M2, 

For the free-particle spinor, u(P, s) — u(P, 3)^74. 

The spinor can be used to construct a positive energy projection operator: 

'^+^^^ '■=^^ '^ = ^ ^-'^ -P + M). (A.7) 

A negative energy spinor satisfies 

v{P, s) (17 • P - M) = = (17 • P - M) v{P, s) , (A.8) 

and possesses properties and satisfies coustraiuts obtained via obvious analogy with u{P,s). 
A charge-conjugated Bethe-Salpeter amplitude is obtained via 

r{k;P) = c'' r{-k;Pf c , (A.9) 

where "T" denotes a transposing of all matrix indices and C = 7274 is the charge conjugation 
matrix, = —C. 

In describing the A resonance we employ a Rarita-Schwinger spinor to unambiguously 
represent a covariant spin-3/2 field. The positive energy spinor is defined by the following 
equations: 

{iy P + M)ui,{P;r) = 0, -fi,Ui,{P;r) = , Pi,u^{P;r) = , (A.IO) 
where r = —3/2, —1/2, 1/2, 3/2. It is normalised: 

u^{P;r')u^{P;r) = 2M, (A.ll) 

and satisfies a completeness relation 

3/2 

— Yl MP;r)u^{P;r)=A+{P)R^,, (A.12) 

r=-3/2 

where 

1 2 - - 1 - 

= S^^Id - g7M7!' + -^PuPi^Iu - i-^lP^-y,^ - P^liA , (A.13) 

with = — 1, which is very useful in simplifying the positive energy Zi's Faddeev equation. 
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